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Abstract 

Conventionally, perturbative and non-perturbative calcula- 
tions are performed independently. In this paper, valleys in the con- 
figuration space in quantum mechanics are investigated as a way to 
treat them in a unified manner. All the known results of the interplay 
of them are reproduced naturally. The prescription for separating the 
non-perturbative contribution from the perturbative is given in terms 
of the analytic continuation of the valley parameter. Our method is 
illustrated on a new series of examples with the asymmetric double- 
well potential. We obtain the non-perturbative part explicitly, which 
leads to the prediction of the large order behavior of the perturbative 
series. We calculate the first 200 perturbative coefficients for a wide 
range of parameters and confirm the agreement with the prediction 
of the valley method. 



In this paper, we report on "valleys" in the configuration space of the path- 
integral. The valley is a series of configurations which connects a classical min- 
imum (vacuum) with an instanton-like configuration. As will be shown in the 
following, it reveals how non-perturbative effects interfere with perturbative ef- 
fects. 

The perturbative expansion is the most common method to analyze models in 
quantum mechanics or quantum field theories. Its predictive power is, however, in 
doubt in a strict sense, for the perturbative series is known to be divergent in most 
of the quantum systems Q. Models in quantum mechanics have been intensively 
studied in order to understand the problem and to obtain physical predictions 
[|J. The symmetric double- well potential provides a very useful testing ground 
in this respect ||. Although the energy spectrum must be real by definition, 
the naive application of the Borel-sum method to the perturbation series results 
in a complex spectrum. The interplay of non-perturbative effect, the tunneling 
between the wells, was inferred to resolve this complexity 0, ||. 

Our analysis sheds a new light on this topic: First, it explains naturally 
why the Borel singularity of the perturbation theory is canceled by the non- 
perturbative contribution induced by instantons, or more generally, valley-instantons. 
The interplay of the perturbative and non-perturbative contribution can be un- 
derstood by the decomposition of the integral contour for the collective coordinate 
of the valley. Second, it predicts the large order behaviors of the perturbative 
series for the asymmetric double-well model, which interestingly defeats widely 
held beliefs that the perturbative series would be Borel-summable for states with- 
out quantum tunneling ||. The valley- instanton [|7|, provides us with the basis 
for the path-integral analysis. Some assumptions that underlie our analysis are 
tested by comparing our theoretical predictions with the numerical and exact 
calculations of the perturbative series. 

The model we study is a one-dimensional system with coordinate q, with the 
following Hamiltonian; 

Id 2 

H = -2i? + v ^ 

V(q) = - gqf - m. (i) 

where g(> 0) is a coupling constant and e(> 0) is a parameter which represents 
the asymmetry of the potential. For eg 2 < V3/18, this potential has two min- 
ima. Although our consideration is not restricted to the special kind of model, 
the reason we adopt the model is that it becomes supersymmetric at the special 
values of e. For e = 1,2, this model is known to be supersymmetric || [LJ], [Hfl . 
The non-renormalization theorem makes clear the relation between the pertur- 
bative contribution and the non-perturbative one and gives a rich structure. We 
will see that the result of the valley analysis agrees with the predictions of the 
super symmetry. 
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In the case of e = 0, the simplest example of the valley represents an instanton- 
anti-instanton pair. It approaches to the infinitely separated instanton-anti- 
instanton configuration at one end and has a vacuum configuration at the other 
end; namely, as the instanton and the anti-instanton approach each other, the 
action decreases smoothly, and finally when they coincide, they annihilate each 
other and the configuration becomes the vacuum. A precise definition of the val- 
ley is given by the valley equation [|12|, [R|. The valley equation which we adopt 
here is 



dr' 



5 2 S 



5S 



A 



5S 



Sq{T)5q{r') 5q{r') 5q(r) ' 



(2) 



where 5* is the action and A is a constant JT3|. The configuration which satisfies 
this equation extremizes the norm of the gradient vector / dr(5S/5q) 2 on the 
contour plane with S fixed. Then, if we identify the value of the action as the 
"height" of the configuration, this configuration is in a "valley" in the configu- 
ration space. The instanton-anti-instanton pair is not a classical solution, but 
it can be defined as a solution of the valley equation ||13|| . When instanton and 
anti-instanton are well-separated from each other, the collective coordinates of 
the valley are the relative distance between them and the center of mass. 

The power of the valley equation is demonstrated in the case of e ^ 0. For this 
case, there is no instanton-like classical solution which starts from one minimum 
of the potential and ends at the other minimum. But such a solution exists for 
the valley equation @) with A set to zero 0. We dub the solution which starts 
from the (right) left minimum as the (anti-)valley-instanton. When eg 2 is small, 
the valley-instanton for the present model may be constructed in a similar way 
as given in M; 



eg + -e^ + ^-Te"+ T 
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9 1 + e~ r 
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- + eg e 

9 9 



eg + - 



3eg _ 
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if r < -1 ; 
if -1/eg 2 < t < 1/eg 2 
+ ■■■ if r > 1 , 



(3) 



where u± = 1 =F 3eg 2 + 



When e goes to zero, the (anti-)valley-instanton 



smoothly converges to the ordinary (anti-)instanton. This gives a natural exten- 
sion of the instanton-anti-instanton pair in the case of e ^ 0. There are two types 
of the valleys, one made of a valley-instanton-anti-valley-instanton pair and the 
other of an anti- valley- instanton-valley-instanton pair 0. We call the former as 
I-A valley and the latter A-I valleyj] 

1 An interesting property of the I-A valley (A-I valley) is that it contains the bounce solution 
if e > (e < 0). In most theories, the bounce solution is a signal for instability, then the 
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Now let us evaluate the transition amplitude from the left minimum |0_l) to 
itself 



lim (0 



-HT 



0l>, 



(4) 



in the background of the I-A valley. For simplicity, we focus mainly on the case 
of e = 0. Except for the collective coordinate of the valley, the path integral can 
be performed by the Gaussian integral. Then we obtain 



Z= lim const, x -J= C dRe- s{R)/9 \ 



(5) 



where R is the collective coordinate of the valley corresponding to the relative 
distance between the valley-instanton and the anti-valley-instanton, and S(R)/g 2 
is the action of the valley.0 A factor Tj^fg 1 results from the integration of the 
collective coordinate for the translational symmetry. When e = 0, the action of 
the valley behaves as ]13| 



S(R) 



d2 

— for R -> ; 
2 

- - 2e~ R for R -»• oo. 
I 3 



(6) 



The valley at R = is the vacuum and at R ~ oo the well-separated instanton- 
anti-instanton pair; the term 1/3 in Eq.(|6]) is their action and —2e~ R is the 
interaction between them. If we change the integration variable R to t = S(R), 
the amplitude © becomes 



lim const, x 



T 



Cv 



dtF{t) e 



-t/9 2 



where x = S(T), Cy = [0,x], and F(t) is the Jacobian. From Eq.( 
that the Jacobian behaves as 



(7) 



it is found 



F(t) 



2t 
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for t 
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1 

3 



(8) 



existence of the bounce solution in our stable system was somewhat mysterious. In our valley, 
the integration of the negative mode of the bounce turns out to be the integration of the 
collective coordinate of the valley which corresponds to the relative distance between valley- 
instanton and anti-valley-instanton and cause no instability in the energy spectrum Q . 

2 When the instanton and anti-instanton are not far enough apart from each other, there is 
no natural definition of R. We define R so as to give Eq.((|) and (jfy . 
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To obtain the full form of the Jacobian, the detailed analysis of the valley is 
needed. Instead, we simply assume that the full form of the Jacobian is given by 



F (t) = sttJ^I IT . (9) 



'2t (1/3 - 1) 

where f(t) is an analytic function that satisfies 

/(o) = ^ /(i/3) = (io) 

No more details of f(t) are needed in the following analysis. 

The integral ([?]) contains both the perturbative contribution at t ~ and 
the non-perturbative one at t ~ 1/3. To separate the perturbative and non- 
perturbative contributions, we deform the contour Cy to the sum of Cp and Cnp 
as is shown in Fig.[TJ. Then the amplitude becomes 

T 



1 f 

Z = lim const, x -= / dtF(t) e~ t/9 

x-*l/3 y/g 2 JCp 



T I' 

+ lim const, x -= / dtF(t) e~ t/9 ' 2 . (11) 



x^l/Z yfg 1 Jc NP 



Note that there is a significant resemblance between the first term on the right- 
hand side and the formal Borel-summation of the perturbation series. The singu- 
larity of F(t) at t = 1/3 is the same as that of the Borel function, which indicates 
the non-Borel-summability of the perturbation series. One of the present authors 
conjectured that this decomposition is essential to understand the interplay of 
the instanton effect to the perturbative calculation |TJ| . We identify that the first 



term of Eq. (fiT|) is the formal Borel-summation of the perturbative series and sec- 
ond term is the non-perturbative one. We show that this decomposition naturally 
reproduce all the known results of the interplay of them. We also present the 
predictions of the conjecture and test them. We denote the first term of Eq.([TT|) 
as Zp and the second term as Z-^p- 

An immediate consequence of the decomposition is 

ImZp + ImZ NP = 0. (12) 

This is because Z = Zp + Znp is real. This simple equation explains why the 
imaginary part of the formal Borel-summation of the perturbation series is can- 
celed by that of the instanton contribution. At the same time, it also shows that 
non-zero imaginary part of the non-perturbative contribution is a necessary and 
sufficient condition for the non-Borel-summability of the perturbative expansion. 
Furthermore, when the imaginary part is not zero, we can predict the large order 
behavior of the perturbative contribution from the dispersion relation. To this 
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end, we examine the analyticity of Z P (g 2 ) in the complex g 2 -plane. When the 
phase of g 2 changes to 2tc, the perturbative amplitude Z P (g 2 ) changes as 



Z P (g 2 ) - Z P (g 2 e 2m ) = Z P (g 2 ) - 2ilmZ P (g 2 ), (13) 

then, Zp(g 2 ) has a cut on the real axis in the complex g 2 plane. This is the only 
singularity near the origin. Thus the dispersion relation becomes 

ry , 2N 1 f°° 1 ImZp(z) 

= -if^f, (14) 
TV Jo z — g 

where we have neglected the contribution from the singularity far from the origin. 
If we formally expand the integrand, we obtain the following;^ 

" rn=0 JU 41 

When m is large enough, the coefficient of g 2m is determined by the singularity 
near the origin. Then, we obtain the large order behavior of Z P (g 2 ) as 

oo 

Z P (g 2 ) = £ z m g 2m , (16) 

m=l 

m ^oo 1 r°° ImZ NP (z) 
^ m ~ / dz —rj—, (17) 

7T JO Z m+1 

which reproduces the formula in |T5|| . 

We find that Z-^p also reproduces the old result for the instanton-anti-instanton 
contribution as follows: For weak coupling, the dominate contribution of the in- 
tegral of .Znp comes from t = 1/3. Then if we expand F(t) as 



and use x = 1/3 — 2e T + • • •, we obtain 

T-ln(-2//)- 7 l, (19) 



T 2 

Z NP = const, x —p^ e~ 1/3g 



where 7 is the Euler's constant. This is the exactly same as the old result An 
advantage of our viewpoint is that it clarifies the meaning of a formal analytic 
continuation in the coupling constant P, H. For this purpose, we perform an 



3 This result holds for a sufficiently large m even if Zp(z) has power divergence for \z\ — > 00 
which necessitates subtractions for the dispersion relation Eq.(|l4|). We would like to thank 
Prof. S. Matsuda for bringing this point to our attention. 
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analytic continuation for Znp(<? 2 ). We denote the contour for Z^p{\g 2 \e l6 ) as 
Cnp(#). If 9 is changed from zero to tt, the contour Cnp(O) changes to C^p(n) 
(see Fig. [J) and in the weak coupling limit, Cnp(vt) can be safely replaced with 
Cy. The resultant integral coincides with that obtained by the formal analytic 
continuation for Z(g 2 ); Z^p(\g 2 \e l7T ) = Z(\g 2 \e l7T ). Therefore, in our viewpoint, 
the formal analytic continuation for Z(g 2 ) is naturally justified as the real analytic 
continuation for Z^p(g 2 ). 

For e ^ case, although some complications arise, the main result does not 
change. Especially, Eq. ([T6|) and ( |17|) also hold in this case. 

In general, there exist valleys which connect the classical minima to the multi- 
valley-instantons. These play important roles in the calculation of the spectrum of 
the excited states. In the following, we incorporate them into the partition func- 
tion, limr^oo Tr(e~ HT ), which we also denote by Z. We have not yet completed 
the analogous analysis of the I-A valley, but the present approach is sufficient to 
evaluate the non-perturbative part Z^p in the weak coupling limit. 

When all the valley-instantons are well separated from each other, the action 
of the multi- valley-instantons is derived by the standard technique ||, |IJ. If the 
valley consists of n pairs of the (anti-)valley-instantons, we obtain 

n 9 n o n 

°y i=i 9 i=i 9 i=\ 

where Ri is the distance between the i-th valley-instanton and the i-th anti-valley- 
instanton and Ri is that between the i-th anti- valley-instanton and the (z + 1)- 
th valley-instanton. When g 2 = \g 2 \e 17T and \g\ 2 <C 1, the following equation 
is meaningful since the interaction between the valley-instanton and the anti- 
valley-instanton is repulsive in this case. (Note that the interaction eJ^Ri is n °t 
attractive since eJ^Ri — e T — e Ri- See Eq. (|22|) below.) 

oo 

Z= lim E« 2 ^n(T), (21) 



n=l 



where a = e 1 ^ g2 /gir 1 ^ 2 , and 

J n{T) = - Y[dR i dR i 5[J2(Ri + Ri)-T) 
n Jo i=i \i=i ) 

^E^ + iE^ + ^E^ • ( 22 ) 

i=i y i=i y i=i / 

One lesson to draw from the analysis of the I-A valley is that Z is equal to the 
non-perturbative amplitude Z^p when g 2 = \g 2 \e %1 * and \g 2 \ 1. Turning back 
to 9 = after the evaluation of Z, we obtain Z^p for positive g 2 . 



6 



We extend the Zinn- Justin method to evaluate Eq.([H|). If we rewrite the 
delta function as 



/ n _ \ 1 rioo-r) ( n „ \ 

5\Y^R l + R l )-Tj = —J^ds e xp^-sT + sJ2(R i + R l )y (23) 



the integrals over Ri and Ri are factorized. After summation over multi-valley 
instanton contributions and turning back to 9 = 0, we finally obtain 

1 r*oo-ij ^d/fs 



np = - lim / dse- Ts -^, (24) 



where 

r°° ~ ( 2 ~2~ N 

= I — a 2 dRdRexp[(s + e)R + —e- R + sR+—e- R 
Jo \ a 1 q 2 



9 9 

f \ 2s+t 

,2 



= l-a^--j r(-s-e)T(-s). (25) 

If we denote the poles corresponding to the solutions of the equation 0(s) = as 
s n , we obtain 



E^ T (26) 



Therefore, s n give the non-perturbative contribution to the energy levels. 

For small coupling, s n can be obtained as a perturbative series in a. For e 
away from integer values, it yields the following to the first nontrivial order, a 2 ; 

E$(e,N) = i + iV + aW(e,iV)a 2 , 

E$(e,N) = -e + i + iV + a(-)(e,iV)a 2 , (27) 



where 



_l)N+i / 2 \ ±£+2Ar 



« (±) ( e ,iV) = ^^r( Te -iV)^-j . (28) 

The energies (j2~7D are for the localized states in the left well and the right well, 
respectively, as can be seen in the free limit, a — > 0. For e = iV"o(= 0, 1, 2, • • •), 
only E^p (e, N) for iV < Nq in the above is valid. The rest has to be solved taking 
into account the confluence of the poles of the two T-functions in Eq. (p5|) , which 
yields the following to the order a 2 ; 



E NP (N ,N,±) = l + N±a 



2 \ 2N+N 



2 \ N\ (N + N )\ \g* ) 

a 2 ( 2 \ 2N+No i / / 2 \ N 1 N+N ° 1 



+ t{7) nhnTW- 21n (-^ +27 -S«-S 



/7 



(29) 
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The plus and the minus signs in the above expression Eq.(|29|) correspond to 
two linear combinations of the perturbative states in the left and right well with 
the same zero-th order energy. This situation is analogous to the lifting of the 
degeneracy by the instanton contribution for the symmetric double-well potential. 

The expression of the nonperturbative contribution to the energy levels, Eq. (|27 
- Eq. (|29"D , contain imaginary parts, which are of order a 2 and are continuous at 
integer e's. Using Eq.(|l7]), these imaginary parts lead to the following leading 
term of the m-th order perturbative coefficient of the iV-th level, E P (e, N,m) 
(E P ±] (e,N) = E^ =0 E P ± \e,N,m)g 2m y, 



E P ±] (e,N,m) = A (±) (e, iV)3"T(± e + 2N + m + 1] 
3 6 ±e+2Ar 



1 + 



A (± \e,N) = 



'vriV!r(±e + 1 + N)' 



, (30) 
(31) 



The 



Ref. 11 



expression 
0- 



for 



E^U, 



1) coincides with the expression obtained in 



We have independently carried out the numerical and exact calculation of the 
perturbative coefficients Ep(e,N,m) by the methods described in Ref. fTT|, p~5 , 
18(1 to the 200-th order for the following four categories; (a) N = (— ) level 
(ground state) for e = to 10 with Ae = 0.2 interval, (b) N = (+) level for 
e = to 20 with Ae = 0.2, (c) iV = 3 (+)-level up to e = 6.5 with Ae = 0.5, (d) 
iV = 1 to 6 (— ) levels for e = 2.5. In order to check the leading m-dependent 
terms in Eq. (|30|) , we take their ratio, 



E^(e, 



N,m) 



E^ ] (e,N,m- 1) 



3(±e + m). 



(32) 



The agreement between the numerical fitting of the perturbative coefficients for 
m = 150 to 200 and this theoretical prediction is excellent in all the calculated 
cases, with the maximum error of order 0.1%. Next we have calculated A^\e, N) 
defined by Eq . (|30"D numerically and compared with Eq. (|3ll) . The result for the 
case (a) is plotted in Fig.[| The difference between Eq.(|3"ID and the calculated 
value is at most 0.1 % (at e = 9.8). For the case (b), the error is at most 15 % (at 
e = 20). For the cases (c) and (d), the maximum error is 0.15 %. In summary, 
the agreement between the expression (|30|), (|3TD, and the actual perturbative 
coefficients is excellent in all the cases examined. 

The reader may note that A^~\e, 0) is zero for any positive integer e. For 
e = 1, this is because of the supersymmetry, which prohibits any perturbative 
correction to the ground state energy 0. For e = 2, there are no perturbative 
corrections to the energy levels of the ground state and the first excited state due 
to the supersymmetry [IT]. We have calculated these energy levels numerically 
and compared them with our valley result, Eq.(|27D - Eq. fl29|) . The result is plotted 



S 



in Fig.|] and the agreement is excellent. For e = 3,4, • • •, there are analogues of 
super symmetry that explain vanishing of and also some other results for 
excited states. These results will be published in near future. 

We stress that the bounce solution plays no important roles in the large order 
behavior of the perturbative series. As a result of this, states that have no as- 
sociated quantum tunneling phenomena do not necessarily have Borel-summable 
perturbative series, in contrast to some beliefs ||. In fact, non Borel-summability 
is found for low-lying stable states when e 7^ integer, which is confirmed by the 
calculation of the perturbative coefficients. 

We believe that our analysis clearly shows that the valley is essential for the 
definition of the so-called "non-perturbative effects"; only when it is defined in 
the light of the separation from the valley, it becomes physically sensible and 
calculable. Same can be said of the perturbation theory; the origin of its Borel 
singularity and its cancellation by the non-perturbative effects become evident in 
view of the valley. 
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Figure 2: The change of the contour C N p(#) as 9 is changed from zero to it. 



11 




Figure 3: Comparison of the fitting to the perturbative coefficients (indicated by 
dots) and the theoretical prediction (solid line) of A^\e, 0). 
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Figure 4: Comparison of our valley results and the numerical results for the 
energy levels of the ground state and the first excited states at e = 2. AE is the 
difference between the energy level and the zero-th order results. The solid line 
is our result, while the dots represent the numerical results. 
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